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Abstract: 



Cn ' Based on macroscopic QED in linear, causal media, we present a consistent theory for the 



Casimir-Polder force acting on an atom positioned near dispersing and absorbing magne- 

^^ ' todielectric bodies. The perturbative result for the van-der-Waals energy is shown to exhibit 

^+ , interesting new features in the presence of magnetodielectric bodies. To go beyond pertur- 

^^ ' bation theory, we start with the center-of-mass equation of motion and derive a dynamical 

expression for the Casimir-Polder force acting on an atom prepared in an arbitrary elec- 

Oj' tronic state. For a non-driven atom in the weak coupling regime, the force as a function of 

time is shown to be a superposition of force components that are related to the electronic 

^ ' density matrix elements at chosen time. These force components depend on the position- 

2 ■ dependent polarizability of the atom that correctly accounts for the body- induced level shifts 

and broadenings. 

PACS: 12.20.-m, 42.50.Vk, 42.50.Nn, 32.70.Jz 



d ■ 1 Introduction 

c^ , 

Being a result of the vacuum fluctuations of the electromagnetic field, Casimir-Polder (CP) 
forces are experienced by any atomic system in the presence of magnetodielectric bodies. 
They play an important role in physical chemistry Tl , and they hold the key to potential ap- 
plications in micro- and nanotechnology such as the construction of atomic- force microscopes 
P] or reflective atom-optical elements [5]- 

On short time scales and for weak atom-field coupling, the CP force is commonly derived 
from the van-der-Waals (vdW) energy calculated by means of time-independent perturbation 
theory 0)13 El or linear response theory [30 EI- Following the former approach, but using a 
quantization scheme for the electromagnetic field in the presence of dispersing and absorbing 
magnetodielectric bodies, we derive a general expression for the vdW energy of an atom 
prepared in an energy eigenstate, where we focus on the influence of the magnetic properties. 



The failure of perturbation theory is evident, because it cannot account for the internal 
dynamics present for an atom initially prepared in an excited state, and because the leading- 
order atomic polarizability, which essentially determines the CP force acting on an atom in 
the ground state, does not incorporate the body-induced energy shifts and broadenings, 
which can become drastic for small atom-surface separations. To overcome this deficiency, 
we present a dynamical treatment by basing the calculations on the well-known Lorentz force 
that governs the atomic center-of-mass equation of motion. The CP force is then obtained 
by taking the average of the Lorentz force with respect to the electronic quantum state of 
the atom and the electromagnetic vacuum, resulting in a dynamical expression that is valid 
for both strong and weak atom-field coupling. For the case of weak atom-field coupling, this 
expression is further evaluated with the aid of the Markov approximation. 

2 Basic equations 

Let us consider a neutral atomic system (e.g. an atom or a molecule) interacting with the 
electromagnetic field in the presence of linear, causal, magnetodielectric media and begin 
with the Hamiltonian 

H = Ha + Hf+Haf- (1) 

According to the multipolar coupling scheme |10j . 

a 

is the Hamiltonian for the atomic system consisting of particles a with charges qa , masses 
rua, positions fa, and canonically conjugated momenta p^, where 



PA(r) = V'fera / dA5(r-fA - Ar-Q 

"'0 



) (3) 



is the atomic polarization relative to the center of mass 



r^^Y. — ^^ (4) 

Tq = f„ - f A (5) 

denoting shifted particle coordinates. The Hamiltonian Hp characterizing the medium- 
assisted electromagnetic field is given by [111 [T^ 

A— e,T?i 

where the bosonic fields fA(r,u;) [and f^(r,u;)], 

/A.(r,c.),/l,,^,(r',c.')] =SxySu'S{r-r')S{cj~u;'), [h^{r,u;), fye{r\u;')\ =0, (7) 

play the role of the dynamical variables of the electromagnetic field plus the medium, where 
A = e and X — m, respectively, refer to the electric and magnetic properties of the medium. 



Finally, the multipolar-coupling Hamiltonian describing the interaction between the atomic 
system and the medium-assisted electromagnetic field in electric dipole approximation reads 



i?. 



AF 



-dE(fA) 



1 



2?TiA 



[pA,dxB(fA)]_^ 



[E(f), electric field; B(f), induction field; [a, b]+ = ab + ba], with 



(8) 



(9) 



being the electric dipole moment of the atomic system. Note that the second term on the 
right-hand side of Eq. (jH)) describes the Rontgen interaction due to the translational motion 
of the center of mass [in| . 

a 

denoting the total momentum of the atomic system. 

The medium-assisted electric field [which in the multipolar coupling scheme has the 
physical meaning of a displacement field with respect to the atomic polarization iPJ] and 
the induction field can be related to the fundamental bosonic fields via ^2 ^] 



E(r) 



B(r) 



duj'E{r, Lu) -FH.c. 



duj (iw)~^ V X E(r, to) + H.c, 



E(r,c.)= Yl fd^''G^{r,r',Lu)h{r',u;). 



The c-number tensors 



Ge(r, r',w) =i' 



ttsq 



Ime(r', oj) G(r, r', uj) 



G,„(r,r',w) = ~i-\ ImK(r',w) G{r,r',uj) x Vr 

C V TTEo 



(11) 
(12) 
(13) 

(14) 
(15) 



where [G(r,r',tt') x Vr'] .. —ejkid'iGik{r,r',uj), K{r,uj) — fi ^{r,uj), are given in terms of the 
(classical) Green tensor, which is defined by the differential equation 



V X K{r,Lu)V X — je{r,LL!) 



G(r,r',tj) = J(r-r') 



(16) 



together with the boundary condition at infinity. Note that the (relative) permittivity e(r, lu) 
and permeability /J,(r, uj) of the (inhomogeneous) medium are complex functions of frequency, 
whose real and imaginary parts satisfy the Kramers-Kronig relations. The Green tensor has 
the following useful properties [TT] . 



G*(r,r',c^) = G(r,r',-a;*), 



(17) 



G(r,r',c.) = GT(r',r,c.), 



(18) 



fd^shiRK{s,Lu) \G{r,s,oj) X Vs [Vs X G*(s,r',w)] 
H — j-Ime(s,cj) G(r, s,Lj)G*(s,r',Lj) ^ = ImG(r,r',a;), 



(19) 



where the last equation can be combined with the definitions H14(l and H15|l to yield the 
useful identity 

J2 J dhGxu'{r,s,u;)GljAr\s,oj) ^^u;HmG,,{r,r',iu). (20) 



A— e,m 



3 The van-der-Waals energy 

Let us consider an atomic system at rest (f a '"^i'a) which is prepared in an energy eigenstate 
\l), i.e., an eigenstate of the Hamiltonian Ha [Eq. |(2J)] written in the form 

HA=Y.E„\n){n\, (21) 

n 

and calculate, within the frame of Schrodinger's perturbation theory, the leading-order en- 
ergy shift 

A^, = A^}°^ + AE^^^ (fa) (22) 

of the state I/)|{0}), where |{0}) denotes the ground state of the fundamental fields i\{r,uj). 
According to Casimir's and Polder's pioneering concept 2], the position-dependent part 
A£'; (fa) of this energy shift can be interpreted as a potential energy 

Ui{rA) = AEl'\rA), (23) 

commonly called vdW energy, from which the CP force acting on the atom in the state \l) 
can be derived according to 

Fi(rA) = -VA;7,(rA) (24) 

CVa = Vr^)- In this approach to the problem, the second term on the right-hand side of 
Eq. (jSl can be disregarded, so that the interaction Hamiltonian that gives rise to the energy 
shift reduces to 

Haf = -dE(fA). (25) 

3.1 General result 

The bilinear form of the atom-field coupling Hamiltonian (|25ll implies that the leading-order 
energy shift is given by the second-order perturbative correction 

^^' = -¥E E ^ / —1- / dV|(/|({0}|-dE(rA)|{lA(r,c.)})|fc)|' (26) 



[V, principal part; ujki = {Ek —Ei)/h; |{lA(r, w)}) =fj^(r,u;)|{0})]. By recalling the definition 
(jll|l together with 113|l . and making use of the commutations relations lO as well as the 
identity H2U|) , it is straightforward exercise to show that 

AEi^^^y^V f duj^^^dikImG{vA,rA,uj)dki (27) 

(d/fc = (Z|d|/c)). In order to extract the position-dependent part of the energy shift in ac- 
cordance with Eq. H22|l . we note that the atomic system should be located in a free-space 
region, where the Green tensor can be decomposed into the (translationally invariant) vac- 
uum Green tensor G^°^ and the scattering Green tensor G^^^ that accounts for the presence 
of magnetodielectric bodies, 

G(r, r', c^) ^ G("' (r, v' ,uj) + G^^^ (r, r', w). (28) 

The vdW energy (|23|l can thus be obtained by making the replacement [G(rA,rA,'^) ^^ 
G(^)(rA,rA,'^)] in Eq. lf77|) . The result can be simphfyed by exploiting the property lfT7|l 
and the well-known asymptotic properties of the Green tensor for large frequencies (cf. 
Ref. jl2| ) to transform the integral along the real frequency axis into an integral along the 
imaginary frequency axis via contour-integral techniques, leading to 

Ui{vA) = Ur{vK) + Ul{vA), (29) 

where 

t/-(rA) = ^ V /°°du^^d;feG(i)(rA,rA,m)dfe; (30) 

is the off-resonant part of the vdW potential, and 

UI{ya) = -^JLo^Q{uJlk)uJfkAlk^eG''''\vA,Yp,,uJlk)AM (31) 

k 

[8(z), unit step function] is the resonant part due to the contribution from the residua at 
the poles at w = iw/fc for loi^ > 0. Introducing the (lowest-order) atomic polarizability 

.,(0)/, X _,. 2^r^ ujki 



ocf^ [uj) = lim - V ^^ d,fe ®dku (32) 

we may rewrite Eq. H3U|) in the more compact form 

Un^A) ^^ r duu^TT[a[°\zu)G^'\rA,rA,iu)]. (33) 

^'^ Jo 

Finally, we note that for an atomic system in a spherically symmetric state we have 

cfHu) = af\u)I = hm A ^ , "^' . |d,,pj (34) 

{I, unit tensor) and thus Eqs. H33() and (|31|l simplify to 

Un^A) = ^ r duu'af'\iu)TTG^'\rA,rA,iu), 



(35) 



C/[(rA) = -^^e(c.ifc)t^ffc|difcpTr [ReG^'\rArA,u;ik)]. (36) 



3 



Equation H29|l together with Eqs. (|31(l and H33(l is an extension of previous results [H] 
to the case of arbitrary causal magnetodielectric bodies. It is worth noting that Eq. H29|l 
also applies to left-handed material |13| . for which standard quantization concepts run into 
difhculties. The calculations presented here can be regarded as the natural foundation for 
similar results obtained on the basis of (semi-classical) linear response theory for dielectric 
bodies 0. 
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Figure 1: The vdW energy Uo{za)12tt^c/{^ouj^ j^d\) of a two-level atom in its ground 
state situated above a semi- infinite magnetodielectric half-space as a function of the distance 
between the body and the interface for different values of ujp^m (wp,e/wT,m = 0.75, ijjT,e/wT,m 

= 1.03, 7e/wT,m=7m/l^T,m = 0.001, Wlo/wT.m = !)• 



3.2 Ground-state atom in front of a magnetodielectric half-space 

To illustrate the influence of the magnetic properties on the CP force acting on an atom near 
a magnetodielectric body, let us apply the theory to a two-level atom [|0), ground state; |1}, 
excited state; cIa =dio =(iA(cos0sin6'e2; -l-sin0sin6'ej^ -|-cos6'e^)] situated above (za > 0) 
a semi- infinite magnetodielectric half-space (2; < 0). Using the appropriate scattering Green 
tensor as given, e.g., in Ref. J^, we find that 



G«(r,r,c.) = 



where 



OTT Jo Po 



fiPo - 



M/3o+/3' 



1 
1 




' p 1 




£/3o-/3 




(37) 
(38) 



denote the reflection coefficients for s- and p-polarized waves, respectively {Pq ^llP' jt? — q^ 
with Im /3o > 0, 0^ = e/iw^/c^ — q^ with Im /? > 0) . Substituting this expression into Eq. H35|l 
and recalling Eq. (|5H) [or equivalently using Eqs. (|S^ and (fM)! and averaging over all possible 
orientations of dA] , we derive the following formula for the ground-state vdW energy (/3o = 



tiaUJiod\ f°° duu'^ f°° ^1 Q „-2bozA 



Uo{za) = ^^^^^^V^ / ^— ^ / ^e— °^- r, 
127r2 7o cjfo + u2 7q 60 1 



a;/ 



(39) 



Based upon a single-resonance permittivity and a single-resonance permeability of Drude- 
Lorentz type, 

2 2 
eiu) = 1 + ^^^^ , ^^{u;) = 1 + ^ ^^^^. , (40) 

Fig.^diplays C/o(za) as function of za for different values of the magnetic plasma frequency 
wp,m- Note that while for short distances the dominant influence of the dielectric proper- 
ties always leads to the familiar attractive 1/zA-potential (cf. Sec. 14. 3|) . strong magnetic 
properties can give rise a repulsive potential barrier at some intermediate distances. 

4 The Casimir-Polder force 

Due to the presence of magnetodielectric bodies the structure of the electromagnetic vac- 
uum can drastically change, resulting in a body-induced shifting and broadening of atomic 
transition lines as the atom comes close the bodies. The perturbative results in Sec. |21 ob- 
viously fail to take such effects into account. Further, the enhanced spontaneous decay of 
an excited atom ^3] leads to dynamical effects that cannot be described within the frame 
of time- independent perturbation theory. In this section, we therefore use a quite distinct, 
non-perturbative approach to the problem, by basing the calculation of the CP force on the 
Lorentz force that governs the atomic center-of-mass motion. 

4.1 General result 

Using the multipolar Hamiltonian (Q together with Eqs. ^, © and ©, and recalling the 
definitions Q and I|1U|) one can easily verify 

WA^A = T[^:"lAfA] =PA + d X B(fA), (41) 

leading to 

mArA = F = ^[i?,PA]+^[dxB(fA)] =|v[dE(r)]+^[dxB(r)]| , (42) 

where in the last step magnetic dipole terms have been dropped in consistency with the 
electric dipole approximation made and terms of the order of 0{v/c) (v, speed of the atomic 
center of mass) have been omitted in accordance with the non-relativistic Hamiltonian (Q . 
Since E(r) and B(r) are defined according to Eqs. 11HI - (|13|I . Eq. 142(1 describes the Lorentz 



force acting on an atomic system in the presence of absorbing and dispersing magnetodi- 
electric bodies and thus generahzes the well-known free-space result |lfij . 

Next, we determine the temporal evolution of the medium-assisted electromagnetic field 
with the aid of Hamiltonian (^ together with Eqs. lHJ, JSJ, and ^. Recalling the definitions 
(|f f|l and 113() and exploiting the commutation relations Q, we may write 

E(r, ^, t) = Ef„„(r, u;, t) + %,^,,,{v, uj, t), (43) 

Ef,3,(r,c.,i)=E(r,w)e— ', (44) 

E,„,_(r,a;,t) = '-^uj' / di' e— (*-*')lmG[r, fA(i'),^]d(i'), (45) 

TT Jo 

where we have used similar approximations as in Eq. H42|) . In accordance with Eqs. (|ll|l - 
(|13|l . we substitute this result into Eq. H42() and take the expectation with respect to the 
field state and the internal (electronic) state of the atomic system, leading to 

(F(t)) = (Ff,ee(t)) + (F,ourcc(i)), (46) 

(Ffrec(t)) = [l "^^ V(d(t)Ef,,,(r, u;, t)) 

+ -:^(d(Ox[VxEf„,(r,w,i)])j +H.C., 



(47) 



i-Ait) 
(Fsourcc(i)) - (Ffourceit)) + {Kol'rc.{t)) , (48) 

X f dt'e~"^'-*-''^V{d{t)lmG[r,rA{t'),u]d{t'))\ +H.c., (49) 

Jo ) r=£-Aft1 



x(d(i) X {VxImG[r,fA(i'),H}d(i'))j +H.c.. (50) 

J r=fA(t) 

Equations H46|I - H5Q(I apply to both driven and non-driven atomic systems and to both weak 
and strong atom- field coupling. In particular, when the medium-assisted electromagnetic 
field is initially in the ground state, then (Ffi.ce(i)) = is valid and Eqs. iB^ and (|Sn|) . 
respectively, just determine the electric part and the magnetic part of the time-dependent 
CP force acting on an atomic system prepared in an arbitrary internal quantum state. 

4.2 Weak-coupling regime 

Representing the electric dipole operator in the form 

d(t) = ^d™„i„„(t) (51) 



{■Amn = \'tn){n\, with |n), \m) being the internal atomic energy eigenstates), we may express 
the dipole-dipole correlation function appearing in Eqs. H49|l and H50|l in terms of correlation 
functions of the Amn (t) as 

(d(i) d{t')) =^Y1 ^™" ® d,r^'n'{A„^n{t)A„^'n'{t'))■ (52) 

m.n ni' ,n' 

In the weak-coupling regime, the Markov approximation can be exploited, and the correlation 
functions {Amn{t)Am'n'{t')) can be calculated by means of the quantum regression theorem 
(see, e.g., Ref. ,17^). For a non-degenerate atomic system with sufficiently slow center-of-mass 
motion this leads to (App. El 

/A (f]A ,,(f'')\-f) ,/A ,('+'^\p{«'^"."(fA)-[r„(fA)-Hr„(fA)]/2}(t-t') /ro^ 

[t>t', m^ n, rj^ ~rj^{t)], where 

Wmn(fA) = UJnin + 5cJm(f a) - (5a;„(fA), (54) 

Ju;™(fA)=^<5wi(fA), (55) 

k 
X k r-^ ^ f^a^ f°°, 2dfemImG(l) (fA,fA,t^)d„jfc 

are the body-induced shifted transition frequencies, and 

r™(fA) = 5]r^(fA), (57) 

k 
r^(fA) = -^e[cJ„ifc(fA)][Wmfe(fA)]^dfe„JmG[fA,fA,'imfe(fA)]d„fe (58) 

are the position-dependent level widths. Note that the position-independent (infinite) Lamb- 
shift terms resulting from G*^°) (f Ai f a, w) [recall Eq. (|28|l ] have been absorbed in the transi- 
tion frequencies uj„in- Equation H56|l can be rewritten by changing to imaginary frequencies 
(cf. Sec. 13.11) . resulting in 

^^ni(^A) = ^e[a)„ife(fA)][w,„fc(fA)]^dfe„ReG(^) [fA,fA,(I'mfe(fA)]d„ife 

Mo /■°°j,_2~ /- ^dfc™G(lHfA,fA,^w)d™fe 



. duu iVkrnirA) TZ .. .„ , — 5 • (59) 

The calculation of 

(i™n(i)) = (Tnrn{t) (60) 

[o-nm(0) = cr„™] yields (App. EI) 

for m ^ n, so the remaining task consists in solving the balance equations 

fT„„„(i) = -r„,(fA)CTmm(i) + ^r™(fA)cr„„(i). (62) 



Upon using Eqs. if^ . l|^ . and ^^, the CP force as defined according to Eqs. (gHJ- 
H5()|l can now be calculated by evaluating the time integrals in the spirit of the Markov 
approximation, resulting in 

(F(t))=^a„™(i)F„„(fA), (63) 

F^„(fA) = F°J,„(fA) + Fr„^(f a), (64) 

with 

™i r \ i^^^X^ f^A 2 V'^d^fcImG(^Hr,fA,cj)dfc„ | 

F™„(rA) = < — > / duju; z — 7:r—. .,-^ ,^ s ;; — ,„ ,-, ,„ > + H.c, (65) 

\^uh cj + cjfc„(rA)-z[rfc(rA) + r„,(rA)]/2 I ^ 

T.mag.~ N [mov- /"""j - ,. . d„feX [VxImG(i)(r,fA,w)]dfc„ 1 

\ T^ , Jq a; + a;fc„(rA) -4rfc(rA) +r,„(rA)J/2 . 

V fe J r— TA 

(66) 

where we have again made the replacement G{r,rx,uj) i— > G^^-'(r, ta,^) [cf- Eq- H28|) ]. As 
f A effectively enters equations (|63f) - (|66|) as a parameter, the caret will be removed in the 
following (ta '^ ta). We finally rewrite Eqs. (|65|l and (|66|) by going over to imaginary 
frequencies (cf. Sec. l3.1|l . Introducing the abbreviating notation 

^mnfe(rA) = Cjnk{vA) + i[r„(rA) + rfe(rA)]/2 (67) 

as well as the generalized atomic polarizability tensor 



Q;m„(rA,w) = -y^ 



h 

k 



Wfe„(rA) - tj - i(rfc[rA) + r„(rA)]/2 

dfen (X) d„ik 



uJkmirA) +UJ + i[rfc(rA) + r„(rA)]/2 



(68) 



with Q:/(rA,w) = Q:/;(rA,w) being the ordinary (Kramers-Kronig-consistent) polarizability 
tensor of an atom in state |/) (cf. Ref. ^Hl)i we derive 

F^^JrA) = F^„"(rA) + F^^^Jr^), (69) 

Fr„*^(rA) = Fr„*^'-(rA) + Fr„^-(rA), (70) 

where 

F^'-(rA) = -[^l duu'[{a^nh{rA,tu) + {a^nh{rA,-tu)]VG[f {r,rA,iu) 

'"""(n) 

KnniT^A) = | A^O JI Q('^"fc)^m„fe (^a) V ® d^feG^l) [r, Ta, r!„„fe (rA)]dfe„ | +H.C., 

(72) 
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Tn'^-{rA)-[\^"fduu^Tr^ 


l^mni'TA) T I ■ \ 


_^™„(rA)^,^(^^^_^^^- 


X [V X G(i)(i",i-A,m)l ) 
/ J 




X (VxG(i)[r,rA,0™„fc(rA)]dfc„)l + H.c. 






J r=rA 



(73) 



(74) 



[(TtT),=Ty,]. 

Equation ^|oi\ together with Eq. (|64|l and Eqs. Ht)9|l - H74() is the natural generalization 
of the perturbative result of Sec. 13.11 In particular, an atom intially prepared in an energy 
eigenstate |^) experiences for Fpt ^ 1 (Fp, characteristic atomic decay rate) the CP force 
(r(i)) ^ (F(0)) =Fl=/(rA) =F;^;-(rA) + F;=;-'-(rA) with 

F^/'°''(i"a) = — -^ \ Auu^\(ai)ij{r K,iu) + (a;),^ (rA, -iu)] VAG'|j^(rA, rA,«u), (75) 
('•A) = ^y]©('^'fe)"?fc(^A)|v®difcG(i)[r,r,r!ifc(rA)]dfc,| + H.c. (76) 



F„ 



[17/fc(rA) = il;;/j(rA); dj/j real]. Note that due to the position dependence of the atomic 
polarizability even the ground-state CP-force FgQ(rA) =FQQ™(rA) cannot be derived from 
a potential in the usual way. It is worth noting that when the atom is initially prepared 
in a coherent superposition of states, then the corresponding off-diagonal force components 
cr„„i (i)Fm„(rA) [n^m) can also contribute to the total CP force as given by Eq. ^^. These 
components contain contributions not only from the electric part of the Lorentz force but 
also from the magnetic part, as can be easily seen from inspection of Eqs. (|73|l and (I74|) . 

4.3 Excited atom in front of a magnetodielectric half-space 

To illustrate the effects of body-induced level shifting and broadening, let us again consider 
a two- level atom in front of a magnetodielectric half-space (cf. Sec. l3.2|l and calculate the CP 
force Fii(2a) acting on the atom in the upper state. For simplicity, we restrict our attention 
to the short-distance limit zx^J\e^\\Lo\/c<^ 1, where the approximations /3o —iq, P — iq can 
be applied to Eq. H37|) . In this limit, which is governed by the dielectric properties of the 
material, Eq. ^^ [together with Eqs. ^^ and ^^] and Eq. ^7^ [together with Eq. (|S5|l ] 
approximate to 



Wio(^a) = t^lO + OWi(za) - OLJq(Za) = WlO ^o,^ fc 3 ur~. ^- M I 112 ' (.' ') 



4(l + COS^g)|£[c^l0(zA)]|'- l 

32TreQhzl |e[wio(2rA)] + 1| 



rfi(l-|-cos^6l) lms[ujio{zA)] , . 

where we have neglected the off-resonant second term in Eq. (|59|l as well as the small free- 
space decay rate. Taking into account that Fii(za) = F^/(za) — F^;''^(za) is valid, from 
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Eq. (jZnj [together with Eq. ^]^ (in the short-distance hmit) and Eqs. IJZZI and IJZHI] we 
finally obtain 



sdKi + cos^e*) \s[nioizA) 



32^0^ \e[nio(zA)] + l\ 
where, according to Eq. H67|l . Hiq^za) =u}io{za) + iTi{zA)/'2. 



1 



(79) 
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Figure 2: The resonant part of the CP force F[;^(zA)AT_gl67r£o x lO'^/^SdA) on a two- 
level atom situated at distance za/At.b =0.0075 of a semi-infinite half-space medium with 
transition dipole moment perpendicular to the interface as a function of bare frequency 
(solid lines), where we have chosen the following parameters; wp^e/wT.e = 0.75, 7e/wT,e = 
0.01, w^ f.dA/i'iTTfieoc^) = 10~^. For comparison, both the perturbative result (dashed fines) 
and the separate effects of level shifting (dotted lines) and level broadening (dash-dotted 
lines) are shown. 

For a single-resonance medium according to Eq. (|40|l . Fig. [3 displays FI^{za) as a func- 
tion of the bare transition frequency wiq. It is seen that the typical dispersion profile of the 
attractive/repulsive CP force around the (surface-plasmon induced) resonance frequency 

-(- ujp g/2, which is already known from perturbation theory, is noticeably re- 



UJS 



- /,.,2 



'T,e 



duced by the decay-induced damping while retaining its width, because the competing effects 
of decay-induced broadening and narrowing due to the frequency shifts almost cancel. 



5 Summary 

Within the frame of electromagnetic-field quantization in the presence of dispersing and 
absorbing linear media, we have presented a consistent theory for calculating the CP force 
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experienced by an atomic system in the presence of an arbitrary arrangement of magnetodi- 
electric bodies. Extending previous results [H] to the case of magnetodielectric media, we 
have derived the vdW energy of an atomic system prepared in an energy eigenstate using 
perturbative methods and demonstrated that magnetic properties of the media can give rise 
to interesting effects such as the formation of a potential barrier for ground-state atoms 
above a magnetodielectric half-space. In an alternative approach based on the Lorentz force 
that governs the atomic center-of mass motion, we have derived an expression for the CP 
force that applies to arbitrary atomic states and both strong and weak atom-field coupling. 
Restricting our attention to a non-driven atom in the weak-coupling regime, we have found 
that the CP force can be written as a superposition of force components weighted by the 
time-dependent intra-atomic density matrix elements that solve the intra-atomic master 
equation. Each force component can be written in terms of the Green tensor for the electro- 
magnetic field and appropriate atomic quantities such as the polarizability. In contrast to 
the perturbative result, the atomic quantities exhibit body-induced shifts and broadenings 
of the atomic transition lines which can noticeably influence the CP force when the atomic 
system comes close to a body. 

Acknowledgement: S.Y.B. acknowledges valuable discussions with O. P. Sushkov as well 
as M.-P. Gorza. This work was supported by the Deutsche Forschungsgemeinschaft and 
the SaxoSmithKline Stiftung. S.Y.B. is grateful for being granted a Thiiringer Landes- 
graduiertenstipendium. 

A Intra-atomic equations of motion 

Upon using the Hamiltonian iQJ together with Eqs. ©, ©, and H21I) . and exploiting the 
decomposition (JSII) the following equations of motion can be derived; 

/•OC 



■ /'OC 

/"OO 

+ / dwEt(fA,w)(d„fcA„A;-dfo„i; 

Jo 





(80) 



where a similar approximation as in Eq. (|42|l has been made. For weak atom-field coupling 
and sufficiently slow center-of-mass motion, the source-field dynamics as given by Eqs. H45|l 
can be approximated by carrying out the time integral in the Markov approximation, re- 
sulting in 

/•CJO 

/ dwE,„^,^„(fA,w) ==^g„„,(fA)imn, (81) 

■^0 m,n 

/•OC 

gm«(rA) = — -I dtJcj^ImG(fA,fA,a;)d,„„C[cj„™(fA)-tj] (82) 

T^ Jo 

[C{x) =ttS{x) +iV/x], with a3„„i(fA) according to Eq. (|5H) . 

We substitute Eqs. H43|) . (|44|l . and (|81|l into Eq. H8U|I and take expectation values with 
respect to the internal atomic state and the vacuum field state, noting that in this case the 
contributions from the free-field part (|44|l vanish. Upon applying the decomposition 



d„fcgfc„(fA) = -iSio'l'^ivA) - irJ;(fA), (83) 
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and using the fact that for a non-degenerate atom ofF-diagonal density matrix elements 
decouple from each other as well as from the diagonal ones, one can derive Eqs. H53|l . (|61|l . 
and linSl. 
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